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ABSTRACT 

Sunyaev-Zcl'dovich (SZ) cluster surveys are considered among the most promising methods for 
probing dark energy up to large redshifts. However, their premise is hinged upon an accurate mass- 
observable relationship, which could be affected by the (rather poorly understood) physics of the 
intracluster gas. In this letter, using a semi-analytic model of the intracluster gas that accommodates 
various theoretical uncertainties, I develop a Fundamental Plane relationship between the observed 
size, thermal energy, and mass of galaxy clusters. In particular, I find that M oc (Ysz/Rsz,2) 3 ^ 4 , 
where M is the mass, Ysz is the total SZ flux or thermal energy, and Rsz,2 is the SZ half-light 
radius of the cluster. I first show that, within this model, using the Fundamental Plane relationship 
reduces the (systematic+random) errors in mass estimates to 14%, from 22% for a simple mass- 
flux relationship. Since measurement of the cluster sizes is an inevitable part of observing the SZ 
clusters, the Fundamental Plane relationship can be used to reduce the error of the cluster mass 
estimates by ~ 34%, improving the accuracy of the resulting cosmological constraints without any 
extra cost. I then argue why our Fundamental Plane is distinctly different from the virial relationship 
that one may naively expect between the cluster parameters. Finally, I argue that while including 
more details of the observed SZ profile cannot significantly improve the accuracy of mass estimates, 
a better understanding of the impact of non-gravitational heating/cooling processes on the outskirts 
of the intracluster medium (apart from external calibrations) might be the best way to reduce these 
errors. 

Subject headings: 



1. INTRODUCTION 

The origin of the present-day acceleration of the Uni- 
verse is arguably the most central question in modern 
cosmology, and is thus likely to dominate theoretical and 
observational efforts in cosmology for decades to come. 
As rece ntly highlighted by t he Dark Energy Task Force 
Report (Albrecht et al. 2006), one of the most promising 
methods for probing the history of cosmic acceleration, or 
its most likely culprit, Dark Energy, is the abundance of 
galaxy clusters at large redshifts, which is exponentially 



sensi t ive to the cosmic expa nsion history (Barbosa et al 
1996 ; Bahcall fc Fan 1998 ). This has motivated many 
upcoming cluster surveys such as APEX, ACT, SPT 
and SZA, which use the thermal Sunvaey-ZePdovich (SZ) 
signa ture ( Sunyaev fc Zcl'dovich 1972; |Carlstrom et al 



[2002 ) of the hot intracluster gas in the microwave sky to 
find clusters at high redshifts 3 . 

However, the accuracy of any cosmological constraint 
inferred from a cluster survey is hinged upon how well 
the mass of the clusters can be estimated from the in- 



dividual cluster observables. For example, Francis et al 
( |2005 ) show that a 10% systematic error in the mass 



estimates is enough to significantly affect the accuracy 
of predicted dark energy constraints from upcoming SZ 
cluster surveys. Although the total SZ flux of a cluster, 
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which traces the total thermal energy of the Intraclus- 
ter Medium (I CM) , is pred i cted t o be a robust tracer of 



its mass (e. 



Si 



irkinshaw 1999: Carlstrom et al. 2005 



Fteid & Spergel p006 ), recent X-ray and SZ observations 
indicate that a significant fraction of cluster baryons may 
have been removed from the ICM, introd ucing a new un 



certai nty into the theo retica l pred i ctions ( Vikhlinin et al. 
" jAfshordi et a.l.j [2005| , |2006| ; |LaRoque~et~ aI| |200C 



2006 



Evrard et al. 2007). Although self-calibration methods, 



through use of phcnomcnological/physical ICM models 



(Majumdar fc Mohr 2004 



ing of clusters ( Lima fc J 1 



Younger et al. 2006), cluster- 
|2004l kOOSlho r gravitational 



lensing ( Bealfon et al. 2006 ; |Hu et al, 2007 ) have been put 
forth as a way to avoid theoretical uncertainties, they do 
rely on ad hoc power-law fitting formulae and/or mod- 
eling assumptions that could jeopardize the accuracy of 
their applications. 

In this letter, I advocate a way to improve the accu- 
racy of mass estimates (or alternatively relax modeling 
assumptions) through including more information about 
the observed SZ profile. In particular, while the usual 
mass estimates only rely on the tota l SZ flux. I develo 



jVClOP 

2002] ) 



a Fundamental Plane relationship ( Verde et al 
among the cluster mass, the total SZ flux, and the SZ 
half-light radius of the cluster. The latter is an indepen- 
dent observable for a moderately resolved cluster, and 
should be readily measurable at similar precisions to the 
SZ flux, for the upcoming SZ cluster surveys. 

2. SEMI-ANALYTIC MODEL OF THE INTRACLUSTER 
MEDIUM 

In order to study the scaling of different ICM observ- 
ables, we first develop a semi-analytic ICM model which 
accommodates a generous allowance for different theo- 
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retical uncertainties. The main ingredient in our semi- 
analytic ICM model is the assumption of hot gas sitting 
in hydrostatic equilibrium in a nearly spherical dark mat- 
ter halo. The dark matter profile is approximated by an 
NFW profile ([Navarro et "aL]|l997|): 



7-200 . We quantify thi s constant through the quantity 
b\ (Evrard et al. 2007), which is defined as: 



b% = 



(T)(<r 200 ) 

M m p(°DM>(< r 200) 



(7) 



p(r) 



Ps 



{r/r s ){l + r/r s y 



(1) 



where r s quantifies the scale at which the slope of the 
density profile changes from —1 to —3. This scale is 
often parameterized using the concentration parameter, 
c 200 = r 200 l r s , where f2oo is the radius within which, the 
mean density of cluster is 200 times the critical density 
of the Universe. We assume a log-normal distribution for 
C200 with the mean: 



where (T) is the mean gas mass weighted tempera- 
ture, /j,m p is the mean particle mass in the ICM plasma 
(/1 ~ 0.59), and (<7q M ) is the mean ID dark matter ve- 
locity dispersi o n. Th e latter is exquisitely constrained in 
Evrard et al. ( [2007 ) through use of a host of different 
dark matter simulations: 



2111! 



0.336 



(C200) 



3.35 



M, 



lQ li h~ 1 M l 



© 



(2) 



and a 22% scatter (Dolag et al. 2004). which is appro- 
priate for an Q m = 0.3 and as — 0-8 cosmology (e.g., see 



K>m) 1/2 (< 7-200) = (1084 km/s, l Mjri , , y 

(8) 

plus 4% random scatter. The 'Santa Barbara Cluster' 
comparison constrains the value of b\ to 0.87 ± 0. 04, for 
a large ra nge of different adiabatic simulations (Frenk 



et al. 199S| ) . More recent high resoluti on simulation s that 



Fig. 11 in Mandclbaum et al. 2006). The NFW gravita- 



include cooling and feedback effects ( Nagai 2006) yield 
consistent values 



tional potential can then be derived analytically 



GAL 



2(ill 



ln(l + r/r s ) 



b% 



0.90 ±0.11, 



(9) 



r ln(l + c 20 o) - c 2 oo/(l + C200)' 



(3) 



Next, we populate this potential with a polytropic gas. 

with T - 1.2 (|5uto et al] 199S|; [AfshOTdi 



et al. 2005; Ostriker et al. 2005). Such a polytropic distri- 



but with larger scatter. We will adopt the latter range 
with a Gaussian distribution. 

In order to set the normalization for the gas pres- 
sure/density, we have to fix the total ICM baryonic 
budget, which we quantif y through f„- a . _ AI^ S (< 



butic m is expected from a turbulent rearrangement and 
i s ro ughly conoioton t with hydro dynamical oimula tion s 



r^no)/ M 2 nn - Various X-ray ( Vikhlinin et alj 2006 \_ Evrard 



et 



( pstriker et al. 200J:) and X-ray observations (e.g., Voit 
et al| 2002 , and references therein). We allow a range 



al. |2007) and SZ observ ations ( |Afshordi et ah 200£ 



2006; |LaRoque et al_2006), as well as hydrodynamical 



i.i< r < 1.3, 



(4) 



with a flat prior, to accommodate uncertainties in and 
deviations from a polytropic distribution. 

In addition to thermal gas pressure, hydrostatic sup- 
port can be p rovided by non-the rmal sources of pressure. 
For example, Nagai et al. (2007) show that subsonic tur- 
bulent pressure can y ield 5% — 20% increase in pressure 



lations ( jNagai " pOOo ) have indicated that / gas may be 
significantly lower than the total cosmic baryonic bu dget, 
which we set to O b /Q m = 0.168 ( |Spergel et aT[|2006| ). To 
accommodate this, we also assume a generous range of: 



0.6 < 



< 0.9, 



(10) 



gradients. Moreover, Pfrommer et al. (2006) argue that 
cosmic rays can contribute up to 32% of the total pres- 
sure in a realistic cluster simulation. To include this 
uncertainty, we consider a wide range of 

5% < 5 nth < 50%, (5) 

with a flat prior, where <5 n th = Pnth/-Pgas, is the ratio of 
non-thermal to thermal pressure components 4 

Plugging all sources of pressure into the equation of hy- 
drostatic equilibrium, and using the polyt ropic relation, 
we fin d the ICM temperature profile (e.g., Ostriker et al 
' 005|): 



for ICM gas mass fractions. 

To estimate the total ICM SZ flux, we need to know 
the outer edge of our ICM model, or the radius of the 
accretion shock, r max . Assuming that gas comes to stop 
at the shock, the temperature behind the shock is roughly 
given by 



1 



MWpV; 2 nf .(l 



(11) 



where Ui n f. is the gas infall velocity 5 . We then use the 
value of infall veloci ty from the spherical collapse model 
QGunn fc GoTtlp^): 



I 



Vint. 



GM{r) 



{V2- [2A c ]- 1 /3 (m/ ^ ) -2/3 



2/3 



T(r) = 



1-r- 1 

1 + 5 nt h 



Mr)-0(r 2 oo)]+T(r 2OO ), (6) 



where T(r2oo) is an integration constant, which is pro- 
portional to the surface pressure of the region within 



0(A, 
( 12 ) 

where A c is the mean overdensity with respect to the 
critical density within the shock radius. Combining Eqs. 
(|6|-(7], [yj|l2|) with mean densities from the NFW profile 
fixes the outer edge of our ICM model. 



4 Note thai. A,u, 'a 
(e.g., pfrommer et al 



200f ; 



expected be 



constant across the ICM 



Nagai et al. 2007). However, the wide 
range ot uncertainty that is already assumed here for 5 nt ^ should 
also include the consequences of its non-uniformity. 



5 Here we have assumed that the non-thermal pressure com- 
ponent behaves as non-relativistic monatomic gas, which is not 
appropriate for cosmic ray pressure. However, we will ignore this 
difference in our model. 
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The final step is to include the ellipticity/triaxiality of 
real haloes in our model. The impact of triaxiality on the 
total SZ flux of a halo is of second order, and so we will 
neglect it in our analysis. However, triaxiality introduces 
a random scatter in the projected SZ profiles, which will 
impact the observed half light radii. To model this, we 
assume that, to first order, the triaxial profile has the 
shape: 



P(r,6,<p) = P(r) 



1 



(13) 



where P(r) is the prediction from our spherical poly- 
tropic model, and Y 2m , s are spherical harmonics. We 
then assume a Gaussian distribution with a reasonable 
range of 

(aL) 1/2 -0.16,( a2m >=0, (14) 

to model the triaxiality of real clusters. This amplitude 
of triaxiality has equivalent moments to an ellipsoidal 
distribution with a xes ratios of 1:0.7:0.5 (expe cted for 
CDM haloes, e.g., Dubinski fc Carlberg 1991), and a 



(spherically averaged) pressure profile P{ 
is roughly con sistent with observat ions and simulations 
of SZ clusters QAfshordi et aT| |2006| ) . 



3. FUNDAMENTAL PLANE OF SZ CLUSTERS 

Let us first quantify the SZ flux of a cluster in terms 
of Ysz, which we define as 

Y sz = Flux(mK • arcmin 2 ) [H (z)d A (z)/c\ 2 hT l (z), 

where Flux(mK • arcmin 2 ) is the total observed cluster 
SZ flux at low frequencies, in units of mK • arcmin 2 , while 
H(z) = IQOh(z) km/s/Mpc and d A (z) are the Hubble 
constant and the angular diameter distance at redshift z. 
Ysz can then be easily described in terms of the proper- 
ties of the ICM model: 



Y sz = 1.022 



M gas h(z) 



(T(keV)). 



(16) 



Now we can generate a random set of 3000 clusters 
uniformly distributed in the range 



13 < log(M 2O o/i/M ) < 16, 



(17) 



with their ICM properties according to the prescription 
that we outlined above 6 . This leads to our mass-SZ flux 
scaling relation: 



200 



1 x 10 i4 M Q //i) Y, 



0.58 

SZ 



(±22% error). (18) 



The error quoted here is the r.m.s. scatter around our 
best fit scaling relation, and reflects a very conservative 
estimate of all the theoretical uncertainties in mass-SZ 
flux relation. Also notice that this includes both system- 
atic and random uncertainties, which cannot be distin- 
guished in our approach. A further simplification is the 
assumed lack of covariance between different uncertain- 
ties. While possible constructive/destructive covariances 
could lead to larger/smaller scatter, their correct account 

6 Notice that none of the assumption that have gone into our 
ICM model would cause a break in the slope(s) of the resulting 
scaling relations, and so the range assumed for cluster masses does 
not change the slope or scatter of the scaling relations. 



4 V* 



• , r .'3 :*.•**• * . 



Eq. shown by s 

mass estimates, Eq. (|19) 
decreases from 22% for 



14 14.2 14.4 14.6 14.8 15 

log (Estimated Mass) 

Fia 1. — Contrast between the usual SZ flux mass estimates, 

ed (black) points and Fundamental Plane 
shown by solid (red) triangle. The error 
former, to 14% for latter mass estimates. 

would require a more detailed understanding of the var- 
ious involved processes. 

To approach the main subject of this letter, i.e. the 
Fundamental Plane of SZ clusters, we will next include 
information about the ICM SZ profile into our scaling 
relation. We do this by calculating the half-light radius, 
Rsz,2, which is defined as the radius of the disk (or cylin- 
der) that contains half of the total SZ flux. The new best 
fit scaling relation for our clusters is: 

A/ 200 = (7.8 x itfHlrJh) Y s °-P (-^|) °' 76 

(±14% error), (19) 

which shows an almost 34% decrease in the error of the 
mass estimate. The new scaling relation is shown by solid 
(red) triangles in Figure (|l|), which should be contrasted 
with the starred (black) points that result from the usual 
mass-SZ flux relation (Eq. |is|). 

One may wonder if using more information about the 
SZ profile of the cluster can help reduce the errors in 
mass estimates even further. To investigate this, we can 
add the radius of the disk containing a quarter of the to- 
tal projected SZ flux, Rsz,i, into the list of observables. 
However, we find that the resulting estimator, which now 
depends on Ysz, Rsza, an d Rsz,a, has a scatter of 13%, 
which is almost the same as the error in the Fundamen- 
tal Plane estimator. Therefore, we conclude that adding 
more details about the SZ profile is unlikely to improve 
the accuracy of mass estimates significantly. 

4. PHYSICAL ORIGIN OF THE FUNDAMENTAL PLANE 

It is interesting to notice that our Fundamental Plane 
relationship (Eq. [l9|) is different from the virial relation, 

Rsz^) 1 ^ 2 



M oc (Y 



SLZ, 



previously adopted e.g., in Verde 



;t alj ( 2002 ). The reason for this apparent discrepancy is 
that the virial relation is only an approximation which re- 
sults from the assumption of hydrostatic equilibrium and 
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TABLE 1 

Breakdown of the fractional error in mass estimates into 
contributions due to different sources of uncertainty. 



sources of AM 2 /M 2 


in Mass-Flux 


in Fundamental Plane 


halo concentration 


0.001 


0.003 


polytropic index 


0.000 


0.000 


gas fraction 


0.004 


0.008 


non-thermal pressure 


0.005 


0.004 


DM velocity dispersion 


0.010 


0.005 


surface pressure 


0.025 


0.009 


triaxiality 


0.000 


0.000 


total 


0.048 


0.022 



self-similar pressure profiles for different clusters. More 
specific information about the initial conditions of the 
cosmological collapse, or the surface pressure, while re- 
laxing the self-similarity assumption, can lead to more 
accurate scaling relations. Since we use the assumption 
of hydrostatic equilibrium, most of our clusters sit close 
to the intersection of the virial relation and Eq. dl9|), 
but (by construction) are better fit by Eq. ([l9|). 

There is a simple way to understan d th e physical ori- 
gin of Eq. (fL9|) analytically. In Sec. 5.1, we show that 



(within our model) the scatter in the mass-flux relation 
is dominated by the uncertainty in the surface pressure, 
which sets the outer boundary of the ICM (Table [l]). 
If we approximate the gravitational potential within the 
ICM region that dominates the SZ flux by an isothermal 
potential, and fix all the cluster/ICM parameters, other 
than its outer boundary, we have M gas oc Rsz,2, while 



(T s . 



const. Therefore, Ysz oc Rsz.2 for a fixed 



cluster potential (and thus fixed M200 ) '■ 

^( A/200) = Ysz I Rsz.i- (20) 
Combining this with the standard scaling relations: 
Y sz oc Mill and Rsza oc M^fo yields: 

Af2oo oc {Y sz /Rszaf /4 - 



(21) 



5. DISCUSSIONS 



5.1. Breakdown of Error Budgets 

In Table we have broken down the errors in mass 
estimates into contributions from different sources of un- 
certainty. To do this for each parameter, we assume that 
it only takes the central value of its assumed distribution, 
and then measure the change in the error quadratures of 
the new scaling relations. 

We see that uncertainties in the ICM gas fraction, non- 
thermal and surface pressures, as well as the dark matter 
velocity dispersion are the main sources of error in our 
mass estimates. In contrast, the uncertainties in the halo 
concentration, the polytropic index, the ICM outer edge, 
or the halo triaxiality have little impact on the errors. 

In particular, the error in both scaling relations seem 
to be dominated by our assumed uncertainty in the ICM 
surface pressure, which is ultimately related with the 
amount of non-gravitational heating/cooling associated 
with galaxy or black hole formation in the clusters. 

Another interesting observation is that the uncertainty 
in gas fraction has a larger contribution to the error in 
the Fundamental Plane relation, than to the error in the 
mass-flux relation. This is due to the fact that Rsz,2 
does not contain any information about / gas (at least 



in our model), while Fundamental Plane masses have a 
steeper dependence on Ysz, and thus are more sensitive 
to /gas uncertainty. 

Notice that possible correlations among different un- 
certainties, that are overlooked in our simple ICM model, 
may tilt the Fundamental Plane from our Eq. (^), and 
also change the size of the errors. However, presence of 
any such correlations is not immediately obvious in our 
current theoretical understanding of the ICM physics. 

One may wonder why the error in our mass-flux rela- 
tionship (Eq. |lj) is so muc h larger than th o se ad vocated 
in numerical studies such as Kravtsov et alj ( [2006 ) , which 
are only ~ 6%. The reason is that these studies measure 
the SZ flux within 3D spheres of fixed ove rdensity (= 500 



times critical, for Kravtsov et al. (B006J)), while our to 



tal SZ flux is integrated out the ICM accretion shock. 
Of course, the latter is a more relevant quantity for 2D 
SZ cluster observations, especially for poor angular res- 
olutions. In fact, our M500 — Ysz,500 relation has only 
a scatter of 9%, which is reasonable as we include more 



theoretical uncertainties than in Kravtsov et al. (2006)'s 
simulations. This shows that the bulk of the scatter in 
our mass-flux relationship comes from the uncertainty in 
the outer edge of our ICM model. 

Finally, we should point out that a breakdown into 
systematic and random errors is also not possible within 
our exercise, due to our poor statistical understanding 
of different non-gravitational processes (such as cosmic 
ray injection or stellar feedback) that affect the scaling 
relations. 



5.2. Noisy Observations 

As the purpose of this letter is to introduce a novel and 
improved mass estimator for SZ cluster surveys, we defer 
a detailed study of the observational issues associated 
with the use of this method to future investigation. Such 
details, while important, should be suited to the specifics 
of each survey, as well as the class of cosmological models 
that one would intend to constrain. However, in what 
follows, I will outline some of the steps that need to be 
taken for a realistic cosmological application. 

We should first recognize that any realistic observa- 
tion of SZ clusters is limited both by the finite detector 
noise, as well as the finite beam resolution. While a poor 
resolution does affect the precision of the SZ flux mea- 
surement, its impact is much more severe for the cluster 
size measurement. For example, if the detector beam is 
significantly larger than the virial radius of the cluster, 
then the Fundamental Plane relation cannot add much 
to the mass-flux relation, even if the cluster is detected 
at several-c level. 

In the absence of perfect resolution, the most practi- 
cal way to use the Fundamental Plane relation is to fit a 
parametrized template (e.g., a Gaussian) to the observed 
cluster SZ map, and replace Rsz,2 with the characteristic 
scale of the template, asz 7 ■ Assuming both Gaussian 
template and beam, this measurement is done by mini- 

7 Of course, the normalization/slopes of the scaling relations 
should be re-calculated for the specific template. Here, we assume 
that the Fundamental Plane or its scatter would not change signif- 
icantly. 
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mizing the following \ 2 function: 



d 2 e 



T 



obs.Z 



Y S ze 



-£ 2 (<r| z+CT g eam )/2 



(27T) 



C t e 



N 



(22) 



where T obs / is the fiat-sky Fourier transform of the clus- 
ter SZ map, Ci is the CMB power spectrum, TV charac- 
terizes the detector noise, and abeam is the size of the 
detector beam. 

In the limit that the detector noise is the primary 
source of measurement uncertainty (C^e _<Tb ° amf <C N), 
the Fisher matrix resulting from Eq. ( p2|) reduces to 
Gaussian integrals which can be calculated analytically. 



In particular, for a well- resolved cluster \o~ sz <; 0b can 



more details about the SZ profile is unlikely to increase 
the accuracy of mass estimates, a better understanding 
of the role of non-gravitational heating/cooling processes 
should significantly reduce the errors. 

While the goal of this letter was to introduce the idea 
of using Fundamental Plane relationship to improve SZ 
cluster mass estimates, there is much more that remains 
to be done in order to exploit the full potential of this 
method. For example, rather than focusing on the cluster 
masses, we can study the cosmological constraints from 
the full bivariate distribution of cluster fluxes and half- 
light radii n(Ysz , Rsz,i)- This combines the traditional 
mass-function constraints, with the constraints res ulting 
from o bserve d scaling relations, as advocated by |Verde| 



the l!ulal (measLimiiuiil+llimiy) enui in Fundamental 

Plane masses is only ~ 28% (22%) smaller than the SZ- 
flux mass estimates for a cluster that is detected at a 
5cr(3cr) level. This is due to the fact that the theoret- 
ical mass degeneracy direction in the Ysz ~ Rsz,2 (or 
Ysz — &sz) plane does not coincide with the degeneracy 
direction of the measured parameters. Including a finite 
resolution will only further deteriorate the performance 
of the Fundamental Plane mass estimates. 

6. CONCLUSIONS 

To summarize, using a semi-analytic model of the intr- 
acluster medium which accommodates different theoret- 
ical uncertainties, we found a Fundamental Plane rela- 
tionship that relates the mass of galaxy clusters to their 
SZ flux and SZ half-light radius. Use of this relationship 
should lead to ~ 34% smaller error in mass estimates 
in comparison to the usual mass-flux relation, and hence 
more accurate cosmological constraints. While including 



et al 



( 2002 ). Another topic that we did not address 
here was the impact of including merging (non-relaxed) 
clusters and/or false detections in our sample. One may 
identify (and thus exclude) these clusters as outliers in 
the Ysz — Rsz,2 plane for a given redshift bin, which 
would not have been possible in the absence of SZ profile 
information. Finally, it is needless to say that the true 
Fundamental Plane, as well as the full impact of differ- 
ent theoretical uncertainties, can only be accurately (and 
adequately) modeled through high-resolution and realis- 
tic cosmological simulations of a fair sample of galaxy 
clusters. 
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